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Abstract The nuclear dynamical deformation, the fusion probability and the evaporation residue (ER) cross sections 


for the synthesis of superheavy nuclei are studied with the di-nuclear system model and the related dynamical 


potential energy surface. The intrinsic energy and the maximum dynamical deformations for “*Cat+*8Cm are 


calculated. The effect of dynamical deformation on the potential energy surface and fusion is investigated. It is found 


that the dynamical deformation influences the potential energy surface and fusion probability significantly. The 


dependence of the fusion probability on the angular momentum is investigated. The ER cross sections for some 


superheavy nuclei in ^Ca induced reactions are calculated and it is found that the theoretical results are in good 


agreement with the experimental results. 
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1 Introduction 


The synthesis of superheavy nuclei (SHN) is an 
important and interesting field in nuclear physics. 
Great achievements have been obtained both in 
experimental and theoretical aspects. In the 
experimental aspect, some experiments with cold 
fusion reactions and hot fusion reactions"?! methods 
were carried out to synthesize the SHN with charge 
number up to 118. Theoretically, in order to 
understand the fusion process and to predict the 
evaporation residue cross sections, several theoretical 
models are put forward, such as the macroscopic 
dynamical modell?, the fusion-by-diffusion (FBD) 


T the nuclear collectivization modell"? 


[15-25 


models! 


the di-nuclear system (DNS) model 
26] 


! and other 
models! The evaporation residue (ER) cross 
sections for SHN decrease with the increase of the 
charge number of SHN, e.g., the ER cross section for 


SHN 118 is less than 1 pb. People are trying to 


produce the SHN with charge number larger than 118 
although in former attempts no expected decay chains 


27,28 
were observed" "??l 


. Thus, it is very important to 
investigate the fusion dynamics and the mechanism. It 
is also of great necessity to give the theoretical 
about the 


combination, the optimal bombarding energy and the 


predictions optimal  projectile-target 
expected ER cross sections of producing SHN because 
only at the optimal conditions could the SHN be 
successfully synthesized for their extremely low ER 
cross sections. If the theoretical predictions can be in 
good agreement with the available experimental data 
systematically, the model may give good predictions 
about the SHN which are not synthesized yet, e.g., 
elements 119 and 120. The heavy ion reaction is very 
complicated where the mass transfer is coupled with 
the dynamical process. Some problems in the fusion 
process are still unknown although much progress has 


been made. For 
[29-31 


example, the dynamical 


deformations l developed in the fusion process 
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should be taken into consideration when one 


investigates the heavy ion reactions. The 
understanding of the dynamical deformation is also 
important in order to calculate the ER cross sections 
for SHN. Thus, paper, 


deformation in the reaction is studied and its effects on 


in this the dynamical 
the potential energy surface and fusion probability will 
be investigated by using the di-nuclear system model 
with dynamical potential energy surface 
(DNS-DynPES model)?! The ER cross sections for 
some "Ca induced reactions leading to SHN with 
charge number 112-118 are calculated and the 
theoretical results are compared with the experimental 
results. The paper is organized as follows. In Section 2, 
an introduction to the theoretical model is given. In 
about the 
fusion probability and evaporation 


Section 3, some results dynamical 
deformation, 
residue cross sections are studied. Finally, a summary 


is given. 
2 DNS-DynPES model 


In Ref.[32], we developed a DNS-DynPES model. In 
this section we briefly mention the formalism for the 
theoretical model. Generally, the production of SHN 
can be separated into three stages: (1) the capture 
process which means the system overcoming the 
Coulomb barrier, (2) the fusion process which means 
the formation of compound nucleus (CN), (3) the 
de-excitation process which means the emission of 
neutrons to take away the excitation energy. The 
evaporation residue cross section in a heavy-ion fusion 
reaction can be written as the summation over all 
partial waves J as 


nA 
Opp (E, ) = 2J+DT(E J 
ER ( TW 2uE,.. »» ) ( c.m. ) (1) 
Fon E on. , J) GR (Eom. 3 J) 
with E.m. being the incident energy in the center of 


mass (c.m.) frame. Here 7(E.m, J) is the transmission 
probability and the survival probability Wou(Ecm., J) is 
calculated using a statistical model. 

The DNS concept was proposed by VolkoyP?! 
and later used to study the competition between 
complete fusion and quasi-fission and to calculate the 


fusion probability in fusion reactions? ^91 


The basic idea of the DNS model is that after 
the capture process, a DNS in the entrance channel is 
formed. Then the DNS evolves via nucleon transfer 
along the mass asymmetry coordinate y instead of the 
direction of the relative distance R between the 
projectile and target. If the nucleon transfer occurs 
from the light nucleus to the heavy one, the mass 
asymmetry 7 increases. Otherwise, if the nucleon 
transfer occurs from the heavy nucleus to the light one, 
the mass asymmetry 7 decreases and the quasi-fission 
rate increases, and that hinders the fusion. Therefore, 
when y is close to 1, a compound nucleus is formed 
with high probability, while the decrease of y away 
from 1, the compound nucleus formation is hindered 
and it may be formed with much lower probability. 

We calculate the formation probability Pen of 
a CN based on the DNS concept. During the nucleon 
transfer process, any configuration of DNS (4), A2) 
with A, =0, 1, ...., Ap + Ar and 45 -Ap*Arz-A, can be 
formed. The evolution of the distribution function of 
each DNS with time can be described by a master 


equation 
aP(A,,t) | 


dt 29s Old, OPA’) 


-d , (t)P(A,0]- A5 OP 2). x 
Since A;+A) -AptAr, only A, is explicitly included 
in the above equation. d,(f) is the microscopic 
dimension for a DNS with a local excitation energy 
Ej, defined in Eq.(3). Ej 
nuclei in this DNS. 

For each nucleus, a valence space is opened 


is shared by the two 


due to the excitation and those nucleons in the states 
within the valence space are active for the transfer 
between the two nuclei. d, (t)= Ca C where Nx 
is the number of valence states and m, is that of 
valence nucleons. AT (f) is the quasi-fission rate of 
the DNS(A;, 42) and W, , (t) -W,, , (t) is the mean 
transition probability between the DNSs(A;, 42) and 
(4 4"). 

In Eq.(2), the microscopic dimension, the 
quasi-fission rate and the mean transition probability 
are all related to the local excitation energy of the 
DNS which reads 


Ejus (4t) =E 


total -— 


Epys(4,t) — Eis (4.0) 3) 
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where 


Es = Eom t (M4 Mj)c (4) 


total 


Ems (4t) = Vins (4t) + (M, +M,)c° (5) 


rot J(J +1) 
E rg je 
a 2I ys (A1) (e 
VyssCA 0.) 9 V,CA 0) +V CA. t). (7) 


Here Vw(41,t) is the nuclear potential calculated with a 
double-folding method?" and Vc(A;,t) is the Coulomb 
interactions potential from the Wong's formula"! In 
this work a tip-tip orientation of the two deformed 
nuclei is assumed. The potential energy in the mass 
asymmetry degree of freedom (the driving potential) at 
t=0, is defined as 


Vos (44) = V, CA f) - Vz (Ast) 


, @) 
+(M,+M,-M,-M,)c’. 


The interaction potential between the two nuclei is 
related to the distance between two centers and in this 
paper the interaction potential takes the minimum 
value in the pocket. 

Because of the attractive nuclear force and the 
repulsive Coulomb force, both nuclei in a DNS are 
distorted and dynamical deformations develop during 


29-31] This results in 


the process of nuclear transfers! 
the time-dependence of the potential energy surface 
(PES). The dynamical deformations of the two nuclei 
satisfy 56°C, / A, = 682C, / A, with the stiffness 
parameter C; (i=1 and 2) calculated from a liquid drop 


modelP?! 


. The total dynamical deformation is the 
average of the dynamical deformations of the two 
nuclei óp = (óP, + 0f, )/ 2. Following Refs.[29, 31] we 
assume that the dynamical deformation evolves in an 


over-damped motion, 
OB, (t) = 8B, (1e "9 ). (9) 


In Eq.(9), the relaxation time m p40x10 ? and the 
maximal dynamical deformation are found when the 


total “intrinsic” energy reaches the minimum, 
Ein (A068) = VAGA By, B.) + VcCA By. B.) 


1 
+> —C.6p?, 
m2 i P, 


(10) 


where the quadrupole deformation /f, = 8° +62, 


(i=1 or 2) with the static deformation parameter p? 
taken from Ref.[40]. 


3 Results and discussion 


The intrinsic energy of the di-nuclear system 
"Ca. Cm as a function of nuclear dynamical 
Fig.l(a) For the 


maximum dynamical 


deformation is shown in 


configuration above, the 
deformation for the system is about 0.32, where the 
minimum value for the intrinsic energy could be found 
and the maximum dynamical deformation is indicated 
as dashed line in the figure. The increase for the 
intrinsic energy of the di-nuclear system when 70.32 
can be attributed to the rises of the nuclear energies 
due to their distortions. In our calculation, it is also 
found that the dynamical deformations increase in the 
beginning 120x10 ? seconds and then gradually get 
saturated. The heavier fragment may get more 
significant dynamical deformation. The interaction 


248 
Cm as a 


potential for the di-nuclear system ^ Car 
function of dynamical deformation is indicated in 
Fig.1(b) and it is found that the interaction potential 


decreases as increasing the dynamical deformation. 
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Fig.1 Intrinsic energy (a) and potential energy (b) as a 
function of the dynamical deformation 6f, for the reaction 
^Ca475 Cm. The maximal dynamical deformation Bmax is 
indicated with the vertical dashed line. 

The potential energy surface (the driving 
potential) for the reaction "Cac Cm is shown in 
Fig.2. The mass and charge distributions evolve from 
the incident channel along the mass asymmetry degree 
of freedom. To form the compound nucleus, the 
di-nuclear configuration must overcome the highest 
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point on the driving potential and the barrier between 
the highest point and the incident channel is defined as 
inner fusion barrier. The solid squares in Fig.2 are for 
the case without dynamical deformation and the open 
circles are for the case when the dynamical 
deformation is fully developed. It is found that the 
potential energy surface changes significantly due to 
the introduction of the dynamical deformation. For 
example, the highest points on the driving potentials 
are 15.4 MeV (for solid squares) and 9.11 MeV (for 
open circles), respectively. The inner fusion barrier 
which hinders the fusion may also get influenced by 
dynamical deformation. The inner fusion barrier for 
the reaction “*Ca+*“*Cm at initial state is 18.93 MeV. 
While the inner fusion barrier is 29.36 MeV when the 
dynamical deformation is fully developed, which 
implies that the fusion probability may be hindered. 


Vees/ MeV 


Fig.2 Potential energy surface (driving potential) for the 
reaction ^Ca4? Cm, 


The fusion probability for the 
^5Ca? Cm at zero angular momentum as a function 


reaction 


of the incident energy is depicted in Fig.3. The solid 
line and dashed line are for the cases with dynamical 
deformation and without dynamical deformation, 
respectively. It is found that the fusion probabilities 
increase with the incident energy for both the cases. 
When the incident energy is below 210 MeV, the 
fusion probabilities increase rapidly while after that 
the fusion probabilities increase slowly. It is also 
found that the 
significantly due to the dynamical deformation. For 
the case of "Ca4" Cm, the fusion probability with 


dynamical deformation is 


fusion probability decreases 


about one order of 
magnitude smaller than that for the case without 
dynamical deformation. In Ref[41], the fusion 
probability for 55Ca? Cm at FEom=210 MeV is found 


to be about 2.5x10 °, which is much larger than the 
result in this work. Part of the difference may come 
from the inclusion of the dynamical deformation in 
this work. And it implies that more investigations 
about the fusion process should be carried out. 
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Fig.3 Fusion probabilities for the reaction ““Ca + Cm at J=0 
as a function of the incident energy in c.m. frame. The solid and 
dashed lines are for with and without dynamical deformations, 
respectively. 

The fusion probabilities for the reaction 
"Ca" Cm as a function of the angular momentum 
(in unit of 7) at three incident energies are depicted 
in Fig.4. The solid, dash dotted and dashed lines are 
for the cases with incident energy being 210 MeV, 
200 MeV and 190 MeV, respectively. The Q value for 
the reaction is about -167 MeV. It is found that the 
fusion probability decreases as increasing the angular 
momentum. In the region with small angular 
momentum, the fusion probability decreases slowly 
while it drops rapidly when the angular momentum is 


larger than about 40 į . 
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Fig.4 Fusion probabilities for the reaction “*Ca+*“8Cm as a 
function of the angular momentum at some incident energies. 
The ER cross sections for some “*Ca induced 
reactions leading to SHN with charge number 
112-118 are calculated by using Eq.(1) The 
experimental optimal ER cross sections, bombarding 


energies and the number of the emitted neutrons ^ *l 
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are listed in Table 1, together with the theoretical 
optimal ER cross sections for the corresponding 
It is found that the 
theoretical results are in good agreement with the 
experimental results. After the SHN 115, the ER cross 
sections for the SHN appear to decrease continuously 


neutron emission channels. 


with the charge number of SHN increasing. 


Table 1 Experimental optimal evaporation residue cross 
sections (in pb) for some reaction channels leading to SHN with 
Z-112-118 and the theoretical optimal results for the 
corresponding neutron emission channels. The incident energies 
in c.m. frame are in the parenthesis (in MeV). 


Reactions Oexp Oih 
8Ca 428 U 2.4577 (194.7,3n) ~ 1.0 (195.0) 
Ca +37 Np 0.9758 (202.9,3n) ~ 2.9 (202.5) 
“Ca 47? Pu 4.4874 (203.7,4n) ~ 2.4 (205.0) 
* Ca 47 Am 8.95! (202.6,3n) ~ 3.5 (200.0) 
* Ca 428 Cm 3.477 (208.87,4n) ~ 4.0 (207.5) 
* Ca +” Bk 1.3775 (211.3,4n) ~ 2.5 (210.0) 
Ca 429 CE 0.571 (210.4,3n) ~ 1.0 (207.5) 
4 Conclusion 
The nuclear dynamical deformation, the fusion 


probability and the evaporation residue (ER) cross 
sections for the synthesis of superheavy nuclei are 
studied using the di-nuclear system model with 
dynamical potential energy surface. The intrinsic 
energy and the maximum dynamical deformations for 
"Cac *Cm are investigated. The effect of the 
dynamical deformation on the potential energy surface 
is found to be significant. The fusion probability as a 
function of the incident energy is calculated and the 
dependence of the fusion probability on the angular 
momentum is estimated. The fusion probability is 
significantly hindered by the nuclear dynamical 
deformations. The ER cross sections for some 
superheavy nuclei in some ^*Ca induced reactions are 
calculated and it is found that the theoretical results 
are in good agreement with the experimental results. 
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